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^ '■ ABSTRACT 

O ; 

t-H ' Aims. The properties of the very high energy (VHE; E>100 GeV) 7-ray emission from the high-frequency peaked BLLac 

C PG 1553+113 are investigated. An attempt is made to measure the currently unknown redshift of this object. 

O ■ Methods. VHE Observations of PG 1553+113 were made with the High Energy Stereoscopic System (HESS) in 2005 and 2006. 

>H+K (1.45— 2.45/^m) spectroscopy of PG 1553+113 was performed in March 2006 with SINFONI, an integral field spectrometer 
of the ESO Very Large Telescope (VLT) in Chile. 
\A , Results. A VHE signal, ~10 standard deviations, is detected by HESS during the 2 years of observations (24.8 hours live time). The 
j_j ■ integral flux above 300 GeV is (4.6 ± 0.6 s tat ± 0.9 sys t) x 10~ 12 cm -2 s _1 , corresponding to ~3.4% of the flux from the Crab Nebula 
C$ ' above the same threshold. The time-averaged energy spectrum is measured from 225 GeV to ~1.3 TeV, and is characterized by 
a very soft power law (photon index of P = 4.5 ± 0.3 s t a t ± 0.1 sys t). No evidence for any flux or spectral variations is found on 
any sampled time scale within the VHE data. The redshift of PG 1553+113 could not be determined. Indeed, even though the 
measured SINFONI spectrum is the most sensitive ever reported for this object at near infrared wavelengths, and the sensitivity 
is comparable to the best spectroscopy at other wavelengths, no absorption or emission lines were found in the H+K spectrum 
presented here. 

Key words. Galaxies: active - BL Lacertae objects: Individual: PG 1553+113 - Gamma rays: observations 



1. Introduction 

Evidence for VHE (>100 GeV) 7-ray emission from the 
active galactic nucleus (AGN) PG 1553+113 was first re- 
ported by the HESS collaboration (jAharonian et al. 2006aP 
based on observations made in 2005. This evidence was 
later confirmed (jAlbert et al. 2007ft with MAGIC obser- 
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vations in 2005 and 2006. Similar to essentially all AGN 
detected at VHE energies, PG 1553+113 is classified as 



a high-frequency peaked BL Lac [(Giommi et al. 1995) and 
is therefore believed to possess the double-humped broad- 
band spectral energy distribution (SED) typical of blazars. 
The low-energy (i.e. from the radio to the X-ray regime) 
portion of the SED of PG 1553+113 is well-studied, in- 
cluding several simultaneous multi-wavelength observation 
campaigns (see, e.g., IQsterman et al. 2006|) . However, the 
only data in the high-energy hump are from HESS and 
MAGIC. The measured VHE spectra are unusually soft 
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(photon index T=4.0±0.6 and r=4.2±0.3, respectively) but 
the errors are large, clearly requiring improved measure- 
ments before detailed interpretation of the complete SED 
is possible. 

Further complicating any SED interpre- 
tation is the absorption of VHE photons 
(jGould fc Schreder 1967] Mauser fc Dwek 2001[) by pair- 
production on the Extragalactic Background Light 
(EBL). This absorption, which is energy dependent 
and increases strongly with redshift, distorts the VHE 
energy spectra observed from distant objects. For a 
given redshift and a given EBL model, the effects of 
the latter on the observed spectrum can be reasonably 
accounted for during SED modeling. Unfortunately, the 
redshift of PG 1553+113 is currently unknowrQ- To date 
no emission or absorption lines have been measured 
from PG 1553+113 despite more than ten observation 
campaigns with optical instruments, including EMMI at 
the NTT and FORS2 w ith the 8-meter VLT telescopes 
( jCarangelo et al. 2003[ ISbarufatti et al. 2006)) . Lower 
limits of z > 0.09 (jSbarufatti et al. 2006|) and z > 0.3 
( Carange lo et al. 2003[ ) were determined from the lack 
of detected absorption/emission lines, implying that the 
effect of the EBL is large in the observed VHE data. The 
absence of absorption and emission lines suggests that the 
non-thermal component of the emission from PG 1553+113 
is largely dominant over that of the host galaxy. This is 
consistent with the fact that, although some hints for a host 
galaxy have been suspected, no clear detection has yet been 
found, even in Hubble Space Telescope (HST) images of 
PG 155 3+113 taken during t he HST survey of 110 BL Lac 



objects (Scarpa et al. 2000) Interestingly, approximately 



80% (88/110) of the BL Lacs initially surveyed by HST 
have known redshifts, of which all 39 with z < 0.25 and 21 
of the 28 with 0.25 < z < 0.6 have their hosts resolved, sug- 
gesting that the redshift of PG 1553+113 is indeed large. 
Recently, u sing a re-analysis of the HST snapshot survey 
of BL Lacs, 
of the abso 



Sbarufatti et al. (2005) claimed the dispersion 
ute magnitude of BLLac hosts is sufficiently 
small that the measurement of the host-galaxy brightness 
allows a reliable estimate of their redshift. With the 
assumption that there is no strong evolution, these authors 
set a possible lower limit of z > 0.78 for PG 1553+113. 
However, adding new STIS -HST data of z > 0.6 B LLacs 
to the snapshot survey, O'Dowd & Urry (2005) found 



on the contrary that host galaxies of BLLacs evolve 
strongly, making any photometric redshift determination 
questionable. A conservative upper limit (z < 0.74) was 
determined (lAharonian et al. 2006a[) from the photon 
spectrum measured by HESS. The same limit was simi- 
larly determined from the MAGIC spectral measurement 
(| Albert et al. 20071) . Using both the MAGIC and HESS 
data, a stronger upper li mit of z < 0.42 was later reported 
(Mazin & Goebel 2007) based on the assumption that 
there is no break in the intrinsic VHE spectrum of the 
object. The range of allowed redshift is clearly large enough 
that the significant effects of EBL absorption cannot be 
reliably removed from the observed VHE spectrum of 
PG 1553+113. This correspondingly makes modeling the 
high-energy portion of its SED unreliable. A clear detec- 
tion of the object's redshift would dramatically improve 



the understanding of PG 1553+113. Further, if it were 
found to be distant, its VHE spectrum could potentially 
provide strong constraints on the poorly-measured EBL 
(assuming a reasonable intrinsic spectrum) and contribute 
to establishing the VHE 7-ray horizon. 

Results from 17.2 hours of new HESS observations of 
PG 1553+113 in 2006 are reported here. In addition, a re- 
analysis of the previously published 2005 HESS data (7.6 
hours), with an improved calibration of the absolute energy 
scale of the detector, is presented. HESS and the Suzaku 
X-ray satellite observed the blazar simultaneously in July 
2006. The HESS results from this epoch are also discussed. 
This will enable, for the first time, future modeling of an 
SED determined from simultaneous observations at VHE 
and lower energies. Finally, results of a March 2006 VLT 
SINFONI spectroscopy campaign to determine the redshift 
of PG 1553+113 are also reported. 

2. HESS Observations and Analysis Technique 

A total of 30.3 hours of HESS observations were taken, typi- 
cally in ^28-minute data segments (runs), on PG 1553+113 
in 2005 and 2006. The pointing of the HESS telescopes 
was offset ±0.5° from the position of PG 1553+113 dur- 
ing the observations. The useful exposure is 24.8 hours live 
time, as 9 of the 66 runs fail the HESS data-quality selec- 
tion criteria. The data are processed using the standard 
HES S calibration (_Aharonian et al. 2004)| and analys is 
tools |(Benbow 2005)| Soft cuts (|Aharonian et al. 2"006aP 
are applied to select candidate 7-ray events. The soft cuts 
are selected instead of the standard cuts since they result 
in a lower energy threshold. Hence they are more appro- 
priate given the very soft VHE spectra previously mea- 
sured by both HESS ( T = 4.0; lAharonian et al. 2006a| 
and MAGIC (T = 4.2; I Albert et al. 2007|) . At the mean 
zenith angle of the observations, 37°, the soft cuts result 
in an average post-analysis energy thresholco of 300 GeV. 
Events falling within a circular region of radius 6 cut =0.1A° 
centered on PG 1553+113 are considered on-source data. 
The background (off-source data) is estimated simultane- 
ously to the on-source data from events recorded within the 
same field of view (F oV) using the Reflec ted- Region metho d 
(Berge et al. 2007)| Equation (17) from|Li fc Ma (1983)1 is 



used to calculate the signifi cance of any excess. Th e energy 



of each event is corrected (Aharonian et al. 2006b) using ef- 
ficiencies determined from simulated and observed muons. 
This correction accounts for long-term variations in the ab- 
solute energy scale due to a changing optical throughput 
of the HESS system. The HESS results have been verified 
using independent calibration and analysis chains. 

3. HESS Results 

The results of the HESS observations for each dark period 
in 2005 and 2006, as well as relevant totals for the observa- 
tions, are given in TableQ] The table shows the dark periods 
in which PG 1553+113 was observed, the MJD of the first 
and last night of observations, the dead time corrected ob- 
servation time, the number of on- and off-source events, the 



z = 0.36 is incorrect. 



Falomo & Treves (1990) demonstrate the catalog redshift of 



The soft cuts were initially named the spectrum cuts in 
Aharonian et al. (2006a). 

3 The threshold has been corrected to account for the de- 
creased optical efficiency of the HESS mirrors. 
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Table 1. Results of the HESS observations of PG 1553+113. 
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Dark 
Period 


MJD 
First 


MJD 
Last 


Time 

w 


On 


Off 


a 


Excess 


Sig 

M 


I(>300 GeV£ 
[10- 12 cm- 2 ^ 1 ] 


Crab 6 

% 


X 2 ,NDF^ 


P(x 2 ) 


04/2005 
08/2005 
04/2006 
07/2006 


53492 
53609 
53849 
53938 


53507 
53614 
53860 
53943 


4.9 
2.7 
7.1 
10.1 


1210 
491 
1811 
2236 


8154 
3462 
12742 
15341 


0.125 
0.125 
0.125 
0.125 


191 
58 
218 
318 


5.5 
2.5 
5.0 
6.7 


4.78 ± 1.18 
10.19 + 5.46 
4.08 ± 1.10 
4.33 ± 0.94 


3.5 
7.5 
3.0 
3.2 


2.7,4 
4.9,2 
5.3,8 
3.1,5 


0.61 
0.09 
0.72 
0.68 


2005 
2006 


53492 
53849 


53614 
53943 


7.6 
17.2 


1701 
4047 


11616 
28083 


0.125 
0.125 


249 
536 


6.0 
8.3 


5.44 + 1.23 
4.22 ± 0.72 


4.0 
3.1 


0.9,1 
0.03 , 1 


0.33 
0.86 


Total 


53492 


53943 


24.8 


5748 


39699 


0.125 


785 


10.2 


4.56 ± 0.62 


3.4 


0.7,1 


0.39 



a The systematic error on the observed integral flux above 300 GeV is 20% and is not shown. 
b The integral flux percentage is calculated relative to the Crab Nebula flux above 300 GeV. 

c The x 2 , de grees of freedom (NDF), and x 2 probability P(x 2 ) are given for a fit of a constant to I(>300 GeV) binned nightly 
within a dark period, or monthly within a year, or yearly within the total. 
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Fig. 1. Distribution of 9 2 for on-source events (points) and 
normalized off-source events (shaded) from observations of 
PG 1553+113. The dashed curve represents the 9 2 distri- 
bution expected for a point source of VHE 7-rays at 40° 
zenith angle with a photon index T = 4.46. The vertical 
line represents the cut on 9 2 applied to the data. 



on/off normalization (a), the measured excess, and the sig- 
nificance of the excess. A significant excess of events from 
the direction of PG 1553+113 is clearly detected in each 
year of HESS data taking. A total of 785 excess events, 
corresponding to a statistical significance of 10.2 standard 
deviations (a), is detected in the complete data set. Figure[T] 
shows the on-source and normalized off-source distributions 
of the square of the angular difference between the recon- 
structed shower position and the source position (9 2 ) for 
all observations. There is a clear point-like excess of on- 
source events at small values of 9 2 , corresponding to the 
observed signal, and the background is approximately flat 
in 9 2 as expected. The peak of a two-dimensional Gaussian 
fit to a sky map of the observed excess is located at a J2000 = 
15 h 55 m 44.7 s ±3.0| tat ± 1.3* yst) 6 320 oo 



ll°ll'41"±53£ at ± 



20" yst . The excess, HESS J1555+111, is coincident with the 



position (aj2o p = 15 h 55 m 43.0 s , (fo ooo = 11°1T24.4") of 
PG 1553+113 |(Beasley et al. 2002)| The upper limit (99% 
confidence level) on the extension of HESS J1555+111 is 
1.1'. 
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Fig. 2. Overall VHE energy spectrum observed from 
PG 1553+113. The dashed line represents the best \ 2 fit °f 
a power law to the obser ved data. The upper limits are at 
the 99% confidence level (Feldman & Cousins 1998) Only 
the statistical errors are shown. 



3.1. Spectral Studies 

The photon spectrum for the entire data set is shown 
in Figure [2] These data can be fit (x 2 of 8.4 for 5 de- 
grees of freedom) by a power-law function dN/dE = 
I (E I 300 GeV)" r ) with a photon index T = 4.46 ± 
0.34 sta t i O.lOgyst. Fits of either a power law with an ex- 
ponential cut-off or a broken power law do not yield sig- 
nificantly better x 2 values. It should be noted that each of 
the three highest energy points, E > 0.6 TeV, in Figure [2] 
have statistical significance less than 2a. However, remov- 
ing these points from the fit does not change T significantly. 

The significant photon statistics from each year of obser- 
vations allows the determination of annual spectra. These 
spectra are shown in Figure [31 The results of fits to the 
data of a power law for both years, as well as for the to- 
tal, are shown in Table fJl The lower and upper energy 
bounds, photon index (r), differential flux normalization at 
300 GeV (I ) , X 2 ! degrees of freedom (NDF) , and x 2 proba- 
bility P(x 2 ) for each fit are given in the table. There are no 
significant changes in the spectral shape during the HESS 
observations. By including one negative point at ^600 GeV, 
and two more positive points at higher energies, the 2006 
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Table 2. Results of the best x 2 fit to the annual and overall spectra of PG 1553+113 measured by HESS. 



Epoch 


E ■ 


E 






r 








Io 






y 2 

A 


NDF 


P(y 2 ) 

L Va. ) 




[TeV] 


[TeV] 












[io- 


11 cm- 2 s 


-l 


TeV" 1 ] 








2005 (AH06)2 


0.185 


0.585 


3.98 


± 


0.63 at at 


± 


0.10 syst 


2.59 


± 0.47 stat 


± 


0.52 sys t 


1.7 


2 


0.42 


2005J1 


0.225 


1.265 


4.01 


± 


0.60 stat 


± 


o.io syst 


5.92 


± 1.19 sta t 


± 


1.18 sys t 


2.1 


5 


0.84 


200~ 


0.225 


0.534 


4.45 


± 


0.48 stat 


± 


0.10 syst 


4.35 


± 0.54 stat 


± 


0.87 sys t 


8.5 


2 


0.014 


Total 


0.225 


1.265 


4.46 


± 


0.34 st at 


± 


0.10 syst 


4.68 


± 0.49 stat 


± 


0.94 syst 


8.4 


5 


0.13 



a The AH06 entry is the previously published HESS result l|Aharonian et al. 2006a)> for 2005. The AH06 entry was not 
corrected for long-term changes in the optical efficiency of the system. 

b The 2005 entry corresponds to exactly the same data as presented in AH06, but with a correction (see text) applied to 
account for optical efficiency changes within the data. The 2006 and total entries also have this correction applied. 
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Fig. 3. Annual VHE energy spectra observed by HESS from 
PG 1553+113. Each line represents the best \ 2 fit of a 
power law to the observed data. Only the statistical errors 
are shown. 
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Fig. 4. Integral flux, I(>300 GeV), measured by HESS from 
PG 1553+113 during each dark period of observations. The 
horizontal line represents the average flux for all the HESS 
observations. For each point the time-averaged V = 4.46 is 
assumed. Only the statistical errors are shown. 



spectrum can be fit over the same range as either the 2005 
or the overall spectrum. This fit yields T w 4.9, explaining 
why the overall photon index is slightly softer than the 2005 
and 2006 indices reported in Table O 



3.2. Integral Flux Studies 

All HESS fluxes (i.e. annual, dark period and nightly 
values) throughout this article are calculated assuming 
the measured time- average photon index of T = 4.46. 
Assuming a different value (i.e. 4.0 < T < 5.0) has less 
than a 5% effect on the flux. The observed integral flux 
above 300 GeV for the entire data set is I(>300 GeV) = 
(4.56 ± 0.62 sta t ± 0.91 syst ) x 10~ 12 cm^s" 1 . This corre- 
sponds to ~3.4% of I(>30 GeV) determined from the 
HESS Crab Nebula spectrurr f(A~haronian et al 2006b) | The 
integral flux, I(>300 GeV), is shown in TableUJfor each year 
of observations, as well as for each dark period. The % 2 and 
corresponding probability shown in the table is for a fit of 
a constant to the data when binned by nights within each 
dark period, by dark periods within a year, and by year 
within the total observations. There are no indications for 
flux variability on any time scale within the HESS data. 
Figures [4] and [5] show the flux measured for each dark pe- 
riod and night, respectively. 



3.3. Effect of the Optical Efficiency Correction 



The data previously published ( Aharoni an et al. 2006a|) for 
HESS observations of PG 1553+113 in 2005 were not cor- 
rected for long-term changes in the optical sensitivity of the 
instrument. Relative to a virgin telescope, the total optical 
throughput was decreased by 29% in 2005 and 33% in 2006. 
These losses are due to the reduced reflectivity of both the 
mirrors and Winston cones, and to the accumulation of dust 
on the optical elements. Comparing (see Table [2]) the spec- 
trum for 2005 determined here to that which was previously 
published for the same data shows significant differences in 
the flux normalization (/ ) but not in the photon index 
(r). Figure [S] illustrates this difference which is a result of 
correcting the energy of individual events for the relative 
optical efficiency of the system, determined from simulated 
and observed muon events, as described earlier. The flux 
measured in 2005 is therefore three times higher than pre- 
viously publisher^. The large difference is due to the steep 
spectrum of the source. For harder spectrum sources the 
effect is much smaller. 



The flux was reported above a threshold of 200 GeV. 
Extrapolating the earlier result, using T — 4.0 as reported, yields 
I(>300 GeV) = (1.78 ± 0.45 stat ±0.35 syst ) x 10~ 12 cm^s" 1 . 
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Fig. 5. Integral flux, I(>300 GeV), measured by HESS from PG 1553+113 during each night of observations. The plots 
a, b, c, and d, represent the April 2005, August 2005, April 2006 and July 2006 dark periods, respectively. The horizontal 
line represents the average flux for all the HESS observations. For each point the time-average T = 4.46 is assumed. Only 
the statistical errors are shown. 
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Fig. 6. Spectrum measured by HESS from PG 1553+113 
in 2005 compared to the previously published version 
(| Aharonian et al. 2006a[) . The integral flux above 300 GeV 
is three times higher. Only the statistical errors are shown. 



3.4. Comparison to MAGIC Results 

Table [3] shows the photon index and integral flux 
for PG 1553+113 measured by both the MAGIC 
(jAlbert et al. 2007[) and HESS collaborations. As the 
data were not simultaneously measured, one might not 



expect the measured quantities to agree due to variability 
arguments. In addition, the spectra from the two instru- 
ments are measured over different energy ranges which 
could result in different spectral slopes if the spectrum were 
curved. Nevertheless, the spectral slope (r) measured by 
both instruments agree within errors in all epochs. Further, 
the flux above 200 GeV measured by both HESS and 
MAGIC in 2005 is consistent^. However the flux observed 
in 2006 by MAGIC is lower than the HESS value. It should 
be noted that in 2006, all the HESS data were taken 
after the final published MAGIC observations, whereas a 
good fraction of the 2005 data from both instruments are 
quasi-simultaneous. 

3.5. HESS Results during Suzaku Observations 



(Mitsuda et al. 2007) 
24, 2006, 14:26 



The Suzaku X-ray satellite 
observed PG 1553+113 from July 24, 2006, 
UTC to July 25, 2006, 19:17 UTC with an av- 
erage efficie ncy of ^53% (Suzaku Observation Log 



|http : / / www .astro.isas.ac.jp/suzaku /index. html . en ) . HESS 
observed PG 1553+113 for 3.1 hours live time on the two 
dates of Suzaku observations. The HESS observation log 



° Previously the 2005 HESS result was ~4 times lower than 
the 2005 MAGIC result. The agreement is solely the result of 
the optical efficiency correction. 



6 F. Aharonian et al.: HESS Observations and VLT Spectroscopy of PG 1553+113 

Table 3. Photon indices and integral fluxes measured by HESS and MAGIC. 



Epoch 


HESS" 

r ~ 


HESS 6 
1(>200 GeV) 
[HT 11 cm^s- 1 ] 


MAGIC 

r 


MAGIC 
I(>200 GeV) 
[10" n cm^s- 1 ] 


2005 
2006 


4.01 ± 0.60 stat 
4.45 ± 0.48 stat 


2.21 ± 0.50 stat 
1.72 ± 0.29 stat 


4.31 ± 0.45 stat 
3.95 ± 0.35 stat 


2.0 ± 0.6 stat 
0.6 ± 0.2 stat 


2005+2006 


4.46 ± 0.34 stat 


1.85 ± 0.25 stat 


4.21 ± 0.25 stat 


1.0 ± 0.4 stat 



a The systematic error on the HESS T is 0.1. 

b The HESS fluxes are extrapolated from the values in Table [1] to above 200 GeV using the time-average P = 4.46. The 
systematic error on the HESS flux is 20%. 



Table 4. HESS observation log during the Suzkau epoch. 4. SINFONI Near-IR Spectroscopy 



Start 


Stop 


Time 


Excess 


Significance 


[MJD] 


[MJD] 


[h] 




M 


53940.74231 


53940.76185 


0.453 


27 


2.5 


53940.76328 


53940.78277 


0.436 


20 


1.9 


53940.78431 


53940.80385 


0.437 


7 


0.7 


53940.80531 


53940.82480 


0.435 


-2 


-0.2 


53941.76336 


53941.78291 


0.444 


27 


3.1 


53941.78439 


53941.80392 


0.444 


14 


1.6 


53941.80542 


53941.82497 


0.443 


7 


0.8 



The determination of the redshift of an AGN is gen- 
erally based upon the detection of emission or absorp- 
tion lines in its spectrum . In the unified model of AGN 



on these dates, henceforth called the Suzaku epoch, is 
shown in Table |4j The log contains the start and stop 
times of the seven HESS observation runs from the two 
dates of Suzaku observations. In addition, the live time, 
observed excess and corresponding significance for each 
run are shown. All runs, except the last, are simultaneous 
with Suzaku observations. The last HESS run begins two 
minutes after the end of the Suzaku pointing. 

During the Suzaku epoch, a total of 659 on-source events 
and 4462 off-source events are measured with an on-off 
normalization of 0.125, corresponding to an excess of 101 
events (3.9cx). Due to the low statistics it is not possible 
to produce a 7-ray spectrum from the data. To provide a 
quasi-simultaneous measurement for future modeling, the 
7-ray spectrum observed during the July 2006 dark period 
is shown in Figure [TJi. The best x 2 fit of a power law to the 
July 2006 data yields r = 5.0 ± 0.7, I = 3.6 ± 0.6 x 10" 13 
cm~ 2 s _1 TeV -1 , and a x 2 of 4.8 for 2 degrees of freedom. 
Only the statistical errors are presented for T and I . The 
systematic errors are 0.1 and 20%, respectively. 

The average flus0 during the Suzaku epoch is I(>300 
GeV) = (5.8 + 1.7 stat ± l-2 sys t) x 10" 12 cm" 2 s" 1 . The flux 
for each of the nights in the Suzaku epoch can be seen in 
Figure [SU. In addition, the flux for each HESS run (~28 
min) during the two nights is shown in Figure [Jj). There 
are no significant variations of the run-wise flux during each 
night, or of the nightly flux during the Suzaku epoch. 



(Urry & Padovani 1995) emission lines are generated by 
ionization from the central source. However in BL Lac ob- 
jects like PG 1553+113 the line-equivalent width of these 
emission lines is reduced, possibly even below detectable 
levels, by the strong beaming of the jet continuum, that is 
caused by the alignment of the relativistic jet with the line 
of sight to Earth. Absorption lines can be produced either 
from spectral features of the stellar population of the host 
galaxy or from intervening halos as in the case of quasars. 
The detectability of stellar spectral features depends in- 
versely on the brightness of the central source. Therefore, 
the probability of detection of any spectral features is likely 
optimal in the near infrared (IR) wavelength range. Here, 
the stellar contribution to the overall continuum emission 
is much larger than at optical or mid-IR wavelengths due 
to the shape of the SEDs of the AGN and the stars. In 
addition, the extinction affecting the stellar light is much 
smaller at near-IR wa velengths than at opti cal (e.g. Ah ~ 
0.17 A y ) wavelengths |(Schlegel et al. 1998) 



In an attempt to detect absorption features from the 
host galaxy or emission lines from the AGN, H+K (1.50- 
2.40^m) spectroscopy of PG 1553+113 was performed 
with SINFONI, an integral field spectrometer mounted at 
Yepun, Unit Telescope 4 of the ESO Very Large Telescope 
in Chile. The observations were performed in seeing-limited 
mode using the 125 mas pixel scale, with a spectral resolu- 
tion of R H+K ~ 1500. The seeing was typically 0'.'7 FWHM 
during the 2.5 total hours of on-source observations per- 
formed on March 9, 2006 and March 15, 2006. 

The ESO D ata Reduction Pipeline Version 1.3 
(Jung et al. 2006) and the QFitsViewLl software are 
used for the data analysis. In the first steps, the images 
are cleaned from bad pixels, flat-fielded and wavelength 
calibrated. A distortion correction is then applied and 
the FoV of SINFONI is reconstructed. The output of this 
procedure is a three-dimensional (3D) cube containing 
about 2000 images of the source in 0.5 nm wavelength 
steps. 

Averaging the cube along the wavelength axis enables 
the effective construction of a broad-band image. H-band 
and K-band images were thus extracted from PG 1553+113 



6 Assuming the poorly-measured F=5.01 from July 2006, in- 
stead of the time-averaged V — 4.46, increases the flux by only 
-3.5%. 



7 Written by Thomas Ott (M PE , Garching) ; 
|http://www.mpe. mpg.de/— ot t / QFits View] 
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(a) Spectrum observed by HESS from PG 1553+113 during observations in July 2006. The dashed line rep- 
the best \ 2 fit of a power law to the observed data. The upper limits are at the 99% confidence level 
(Feldman & Cousins 1998)) (b) Integral flux, I(>300 GeV), measured by HESS from PG 1553+113 in each run taken 
during simultaneous Suzaku observations. The horizontal line represents the average flux for all the HESS observations. 
For each point in the light curve the time-average T = 4.46 is assumed. For both plots, only the statistical errors are 
shown. 



data. These images are spatially unresolved and no under- 
lying host galaxy is detected. 

A spectrum is extracted from the final 3D cube by 
summing the individual spectra within an aperture of 5- 
pixel radius centered on the brightest peak in the im- 
age. Observations of two G2V standard stars (HIP071527, 
HIP088235), near both in time and air-mass to the target 
exposures, are used to correct for strong atmospheric (tel- 
luric) absorption. By dividing the science target spectrum 
by the standard star spectrum, the atmospheric absorp- 
tion effects are minimized (see, e.g., IMaiolirio c t al. 19961 
IVacca et al. 2003|) . However, this step introduces features 
into the resulting spectrum that are intrinsic to the 
standard star. Here, these features are accounted for by 
multiplication of the normalized high-spectral-resolution 
atmospheric-transmission-corrected solar spectrurr@ that 
has been degraded to the SINFONI spectral resolution. 

Flux calibration at infrared wavelengths is somewhat 
more problematic than at optical wavelengths since no stan- 
dard spectral catalog exists. A theoretical continuum spec- 
trum of a G2V star could in principle be used to correct 
for the throughput of the spectrograph. However, the flux 
distribution of such stars, as well as later-type ones, in 
the infrared is not well-reproduced by theoretical model- 
ing, and a small error in the stellar parameters (e.g., the 
effective temperature or the metallicity) can yield a large 
error in the continuum. Instead the wavelength-dependent 
flux calibration is determined from observations of two in- 
trinsically featureless (apart from strong hydrogen lines) 
B stars (HIP088201, HIP085812). Here the hydrogen lines 
are manually fit and removed from the B-star spectrum, 
and the result is then divided by a blackbody spectrum at 
the effective temperature of the stars (18500 K and 11900 
K, respectively) to correct for the continuum slope of each 
standard star. As the H and K bands are in the Rayleigh- 



8 NSO/Kitt Peak FTS data of the sun used here were pro- 
duced by NSF/NOAO. 



Jeans regime of the blackbody spectrum, the steepness of 
the continuum slope for these B stars is not particularly 
sensitive to small errors in the effective temperature. The 
final step of the flux calibration is setting the absolute flux 
scale. IR aperture photometry of PG 1553+113 was already 
performed but the source is known to be variable. Therefore 
2MASS H magnitudes of the standard stars, converted to 
flux density, are used instead to scale the total flux observed 
here. This value, 12.7+ 1.3 mJy at 1.65/xm, is app roximately 
half of that observed by Bersanelli et al. (1992)|and is co m- 
parable with R-band fluxes observed (jAlbert et al. 2007)) in 
March 2006. It should be noted that the latter (March 15, 
2006) of the near-IR observations reported here were per- 
formed two days before an outburst in the R-band reported 
bv I Albert et al. 20071 

The final H+K-band spectrum of PG 1553+113 is 
shown in Figure [5J The signal-to-noise ratio is ~250 in 
the H-band and ~70 in the K-band. These measurements 
are the most sensitive ever recorded in the H- and ri- 
bands, and are comparable to the most sensitive spec- 
troscopy of PG 1553+113 at other waveleng ths (see, e.g, 
ICarangelo et al. 2003) ISbarufatti et al. 20061) . The lack of 
data between 80— 1.95 /jm is a result of highly-reduced 
atmospheric transmission between the H and K bands. The 
observed near-IR spectrum is featureless apart from some 
residuals from the atmospheric correction. Thus, in neither 
the broad-band images nor in the spectrum are the influ- 
ences of the gas of a host galaxy or the AGN detected, even 
though PG 1553+113 is bright in the IR. As a result, a red- 
shift determination from these observations is not possible. 



5. Discussion 

The presence of Suzaku X-ray observations, simultaneous 
with the HESS measurements, greatly improves the pos- 
sibilities for accurately modeling the underlying physics 
in PG 1553+113. However, the observed VHE spectrum is 
very soft (r = 4.46 ± 0.34), which could be in large part 
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Fig. 8. Top: H+K-band spectrum of PG 1553+113 extracted from a 5-pixel radius aperture. The gap is due to the highly 
reduced atmospheric transmission between H and K bands. Bottom: Telluric spectrum extracted from the standard stars. 



due to the redshift-dependent absorption of VHE photons 
on the EBL. Modeling the VHE portion of the SED of 
PG 1553+113 requires the effects of this absorption to be 
removed. Unfortunately the redshift of the object is still 
unknown, despite the very sensitive new measurements pre- 
sented in this article. Therefore, modeling of the VHE emis- 
sion from PG 1553+113 is extremely difficult given the large 
range of possible intrinsic spectra. 

PG 1553+113 is not the only bright BLLac for which 
the underlying host galaxy is not detected, even though 
deep, high signal-to-noise ratio observations were con- 
ducted (e.g., O'Dowd & Urry 2005. ISbarufatti et al. 2006]) 
for many. The properties of their host galaxy or their AGN 
may be such that it is presently not possible to derive their 
intrinsic nature and hence redshift. In these cases, it is pos- 
sible that the properties of the host are not those of gi- 
ant elliptical galaxies, but correspond to the faint end of 
the elliptical luminosity distribution (M« ~ —20) or that 
the AGN may even be hosted by a dwarf galaxy, making 
these BLLac objects the radio-l oud equivalent to th e low- 
mass black-hole Seyfert galaxies (Barth et al. 2005) If this 
is the case, it will be very difficult to derive a firm redshift 
for PG 1553+113. 

As mentioned in the introduction, it is possible to derive 
upper limits on the object's redshift using VHE measure- 
ments, some basic assumptions about the source's intrinsic 
spectrum, and an EBL density model. As yet, none of the 
limits derived in this fashion are particularly constraining. 
The strongest, and also based on the most assumptions, 



is z < 0.42 |(Mazin fc Goebel 2007)| Following exactly the 
same methodology as in Aharonian et al. (2006a), the VHE 
spectrum measured here from PG 1553+113 limits the red- 
shift to z < 0.69. 



6. Conclusion 

With a data set that is ~3 times larger than previously 
published ([Aharonian et al. 2006a) . the HESS signal from 
PG 1553+113 is now highly significant (~10cr). Thus, the 
evidence for VHE emission previously reported is clearly 
verified. However, the flux measured in 2005 is now ~3 
times higher than first reported due to an improved cal- 
ibration of the absolute energy scale of HESS. The sta- 
tistical error on the VHE photon index is now reduced 
from ~0.6 to ~0.3. Nevertheless, the error of 0.34 is still 
rather large, primarily due to the extreme softness of the 
observed spectrum (r = 4.46). The total HESS exposure on 
PG 1553+113 is ~25 hours. Barring a flaring episode, not 
yet seen in two years of observations, a considerably larger 
total exposure (~100 hours) would be required to signifi- 
cantly improve the spectral measurement. This large expo- 
sure is unlikely to be quickly achieved. However, the VHE 
flux from other AGN is known to vary dramatically and 
even a factor of a few would reduce the observation require- 
ment considerably. Should such a VHE flare occur, not only 
will the error on the measured VHE spectrum be smaller, 
but the measured photon i ndex may also be harder (see, 
e.g., lAharonian et al.~2 002). Both effects would dramati- 
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cally improve the redshift constraints and correspondingly 
the accuracy of the source modeling. Therefore, the VHE 
flux from PG 1553+113 will continue to be monitored by 
HESS. In addition, the soft VHE spectrum m akes it an ideal 
target for the lower-threshold HESS Phase-Il [(Punch 2005) 
which should make its first observations in 2009. 
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